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THERMAL STRESS IN LONG CYLINDRICAL SHELLS DUE TO 
TEMPERATURE VARIATION ROUND THE CIRCUMFERENCE, 
AND THROUGH THE WALL! 


By J. N. GoopiErR? 


Abstract 


The thermal stress in thin-walled cylinders of any cross section has been 
investigated for internal and external temperatures each varying in any manner 
round the circumference but not in the axial direction. The thickness also may 
vary round the circumference. 

A method is given for calculating the stress from given temperature dis- 
tributions, whatever the shape of the cross section. The stress is evaluated for 
uniform, but different, inside and outside temperatures. 

The circular cylinder is treated in detail and the stress found for the general 
case of circumferential variation. It is shown that the maximum stress will 
depend only on the temperature distributions and the material, and not on the 
thickness or diameter of the cylinder. 


1. The problems of thermal stress, that is, stress due to non-uniform heating, 
like the ordinary problems of stress and strain in elastic bodies, fall into two 
groups. The first group can be satisfactorily treated only by solution of 
the fundamental differential equations of elasticity for appropriate boundary 
conditions.* The second group consists of those problems which can, for 
practical purposes, be adequately treated by the approximate but simple 
and convenient formulas of the subject usually described as the strength of 
materials. The well-known Bernoulli-Euler formula for the bending of 
beams is typical of these. The simplification arises from the circumstance 
of one or two dimensions of the shapes considered—beams, thin plates, and 
thin shells—being small. 


Some of these simple formulas will be employed to solve the problem 
described by the title of this paper. Temperatures varying along the axis, 
but not round the circumference, and temperatures uniform, but different, 
inside and outside a cylinder, have been dealt with by previous writers (4, 6), 
but only for the particular case of the circular tube. The thin tube of any 
shape, with completely general variation of temperature over the surfaces, 
presents a problem of great analytical difficulty. The differential equations 
involved are not linear. The same is true of the flat plate with general 


1 Manuscript received December 30, 1936. ; ] 
Contribution from the Department of Engineering and Metallurgy, Ontario Research 


Foundation, Toronto, Canada. 
2 Research Fellow in Applied Mechanics, Ontario Research Foundation, 


* This group is treated in a forthcoming paper ‘‘On the integration of the thermo-elastic equations.” 
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surface temperatures. The restricted problem we are concerned with, how- 
ever, is not only important in engineering, notably in its application to boiler 
tubes, but also tractable mathematically. 

We begin with the incomplete cylinder, i.e., the section AQPFC, Fig. 1, 
is not a closed curve. The results obtained are necessary for the com- 
plete cylinder, which, for the present, is supposed converted into an incom- 
plete one by a cut along a generator, equivalent to a gap of zero width. The 





Fic. 1. 


problem is reduced, by a physical argument, to the solution of a set of 
equations. A narrow strip PQSR is considered in isolation. It is given 
arbitrary inside and outside temperatures, T; and T, respectively (the 
transition through the thickness being linear), while the remainder of the 
tube is kept unheated. The strip expands, bends and curls, and no longer 
fits the place from which it was taken. But by the imposition of suitable 
terminal forces and couples it is so far restored to its original stete that it 
can be supposed reattached to the rest of the tube, so that points originally 
contiguous are so again without any necessity for straining the rest of the 
tube. The reassembly, however, involves relative rigid body displacements 
of the two cold parts PFC and QA, and hence an alteration of the gap AC. 
This alteration is calculated. Then the stress system required to re-close the 
gap AC, if the tube be a closed one, is found, and the effects of removing the 
imposed terminal forces and couples of the strip are considered. The effects 
of heating over more than an elementary strip follow by integration. 


The Incomplete Cylinder ; 

2. For the moment, let the strip PQSR (Figs. 1 and 2) of (the complete) 
axial length / and infinitesimal arc length 6s, have the inside temperature 
T;, and outside temperature 7., the rest of the cylinder being kept at tem- 
perature zero. 

Let this strip be freed from the rest of the cylinder by cuts along PR and 
QS. Then it is free to expand, without thermal stress.* It increases its 
linear dimensions, and also changes its curvature on account of the tem- 


* There is a general theorem that temperature distributions linear in cartesian co-ordinates give 
rise to no stress (8, p. 204). 


as 
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perature gradient through the thickness. 
If the mean temperature 7, = 3(7; + he 
T.), then the length increases by a uniform ea 
linear strain a7,,, 1.e., by laT, , a being 
the coefficient of linear expansion. The 
arc length 6s increases by 6s. aT». 
The thickness # increases by haT,. 
The length /, which originally had no 
curvature, becomes curved, as sketched 
in Fig. 2 (with JT, > 7;), to a radius r, where 1/r = (T. — Tya/h. 
The arc 6s changes its curvature by the same amount. 


Cs) 





Q Fre. 2. 


3. A compressive stress EaT,, , applied to the ends of the strip will exactly 
cancel the thermal strain aT,, in the direction / (where E is Young’s modulus). 
It will increase 6s by a strain gaT,, (where ¢ is Poisson’s ratio) so that the 
total increase of 5s becomes 6s(1 + a)aTm. It will also increase the thick- 
ness h by hoaT». 

A bending moment of suitable magnitude, applied to the ends, will exactly 
cancel the thermal curvature of /, (T. — T;)a/h. It is easily found from simple 
béam formulas that the necessary bending moment is ¥zEah?(T,—T;)6s. The 
vector representing this moment is in the direction of 6s. 

This moment will produce an anticlastic curvature of 6s, of amount 
o(T. — T;)a/h, and the total change of curvature of 6s is therefore now 
(1+ o)(T. — Ti)a/h. 

The edges PR and QS of the strip are at present thicker than the thickness 
h of the cold tube by h(1 + o)aT,,. In order to bring the faces of the cuts 
of the strip and of the rest of the tube to 
the same configuration, it must be supposed 
that shearing and normal tractions as indi- 
cated in Fig. 3, and :erhaps others, act on 
the strip, and equal and opposite tractions 
on the rest of the tube. 

It is however neither practicable nor 
necessary to compute these tractions and Fic. 3. 
their effects. The strip, at present, is 
regarded as hot, and the rest of the tube as cold. But later the effects of 
heating all such strips will be added so as to arrive at a tube heated all round, 
the temperature varying smoothly. Evidently the adjustment of the thickness 
of one strip to that of its neighbors involves only the temperature differentials. 
In the limit, therefore, the required tractions on the cuts such as PR and QS 
will vanish. We may say that they disappear in the process of integration.* 





* The validity of this argument may be better seen by considering a problem of plane stress. 
Let a long strip be compressed on finite opposite segments on the two edges. The part between these 
segments may be considered in isolation, and we may inquire what tractions are necessary to fit 
it back into the strip. We should find that an adjustment of thickness is required, just as in the 
present problem. But if the compressive forces have a smooth distribution over a succession of 
infinitesimal segments, the process of integration involved eliminates the tractions that have to be 
introduced to adjust the thickness, and the stress is plane. This is confirmed by the exact analytical 
solution of the plane stress problem. 








CANADIAN JOURNAL OF RESEARCH. VOL. 15, SEC. A. 


on 
Nw 


4. For the present, suppose that the tractions act, so that the thicknesses 
become equal. Any other possible effects due to the tractions on the strip 
need not be considered. Then we may suppose that the strip is rejoined to 
the other two parts of the cylinder. But on account of the alterations in the 
length and curvature of 6s, these parts will no longer have their former 
relative positions. 

Consider any cross section of the cylinder and draw rectangular axes 
Cx, Cy with origin on CD. In Fig. 1 these are drawn on the end AQPFC, 
for convenience. 

Then the displacement of AQ relative to CFP has x and y components 


given by, 
~(1 + c)aTn cos. bs — (1+ a)a(T. — Ti) = bs 


and 


~(1 + o)aTm sind. ds + (1 + o)a(T. — Ti) = bs 


The first terms arise from the extension of 6s, the second terms from the 
change in its curvature; the rotation (anticlockwise) of AQ relative to CFP 
is —(1 + o)a(T. — T,)h~és; x and y are the co-ordinates of 6s, and } 
is the angle 6s makes with the x-axis. 

5. If now all the other strips into which the cylinder may be divided are 
treated in the same way, each being given its proper temperatures, the dis- 
placements 6, and 6, and rotation ‘w of the edge AB relative to the edge CD 
may be obtained by superposition. Then 


¥ 
545: 


6, = —(1 + a)af Tx cos $. ds — (1 + oa (7. - Ti) = 


x 


6, = -(it+ oa | T sing@.ds+ (1+ aya (r, — T;) “ds, (1) 


o = -( + o)af (2. = Tr), 


Tm, Te and T; being now given functions of s, the arc length. The integrals 
are taken all round the section. 

These displacements and rotation are maintained by the temperature 
distribution together with: 

(i) Stress —EaT on the ends of the cylinder; 

(ii) A distribution of moment on the ends, 

m = txEah?(T. — T3), 

per unit length of arc, the vector of m being directed along 6s or along the 
tangent at any point of the end section. 

The state of stress anywhere in the wall of the cylinder is given by the same 
stress and moment. There is no stress other than the axial component. 
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6. The end stresses (i) and (ii) of Art. 5 may be eliminated by superposing 
the ordinary isothermal stress distribution having equal and opposite end 
stress and moment. The complete determination of this distribution is a 
problem in the theory of thin shells. However, if the cylinder is long, the 
distribution in the middle part will depend, according to the Principle of 
Saint-Venant, only on the resultant force and couple on the ends. The 
resultant force is EafhT,,ds and the resultant couple has components 


M, = —vzEafh?(T. — T;) cos¢. ds, 
M, = —1zEafh?(T. — T;) sing. ds, (2) 


and from these the corresponding stress in the middle part can be found from 
the simple beam theory. 


The Complete Cylinder 

7. If a cut is made along a generator, the formulas (1) give the displace- 
ment and rotation of one face of the cut relative to the other, when the end 
actions (i) and (ii) of Art. 5 are applied, and the cylinder is heated. The 
integrals are now taken around the circumference. 


The faces of the cut may be brought together by a suitable force and couple 
on one face, an equal and opposite force and couple on the other face, and 
rejoined. This introduces an additional stress distribution to be superposed 
on that of Art. 5. 


To evaluate it, consider the forces X, and Y, and moment M,, per unit 
axial length, distributed uniformly along the 
faces of the cut, acting on the cut tube 
(Fig. 4). The displacements and rotation of 
the face A relative to the face C are con- 
veniently found by Castigliano’s theorem (7, 
p. 434). Since we are dealing with a cylinder 
and not a ring, we use the plate modulus* 
E/(1 — o?) instead of Young’s modulus, £, 
and write D for Eh®/12(1 — o?). The thick- 
ness / will now be taken as constant, although 
there would be no difficulty in retaining it as 
avariable. Then, 


Dé, = M.fyds + X.fy'ds — Yi fxyds } 

Dé, = —M.fxds — X,fxyds + Y.fx*ds (3) 

Dw = M.fds + X.fyds — Y.fxds 
The integrals can be expressed in terms of the mean circumference of the 
section, L, the co-ordinates of the centroid ¥, y, the moments of inertia J,,, 
I,,, and the product of inertia J,,, of the area of the section, since L = fds, 
Ex = /xds, Ly = J yds, ang). = hf yds, yy = hf x*ds, I, = hf xyds. The 
displacements and rotation given by Equations (3) must be equal and opposite 





* Expressing the effect of preventing anticlastic curvature. 
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to the displacements and rotation given by Equations (1). The forces 
X,and Y, and the couple M, can be found by solving the equations, 
M, cage x. ‘ol es — YJz = D(ite)alh{Tncos od. ds+f(Te—T;)yds] , 
— MLhE-—XTry+Volyy = DA+o)alhf[Tn sin d. ds—f(T.—T;)xds], }(4) 
(M, + me — Y,x)Lh D(ito)af(T.—T;)ds . 


ll 


8. The axial stress is made up of four contributions: 
(i) Astress —EaT,, (Art. 3); 


(ii) A stress corresponding to the distribution of an axial bending moment 
tz Eah?(T. — T;) per unit arc length on the mean circumference (Art. 3); 


(iii) The axial stress due to the distribution of traction EaT,, on the ends. 
This cancels, on the ends only, the stress (i) ; 


(iv) The axial stress due to the distribution of the couple on the ends 
—yrEah?(T. — T;) per unit arc length on the mean circumference. This 
cancels, on the ends only, the stress (iii). 


The ends are then completely free. 


The resultant moments of (iv) are given by Equations (2), and the stress 
due to them, in the middle part, 7.e., away from the ends, can be calculated 
by the simple beam formulas. The appropriate moments of inertia are 
not J,, and J,, as defined previously, since the latter are referred to axes 
not through the centroid of the section. Writing I’,,, I’,, for the centroidal 
moments of inertia, and ¢ for the distance normal to the middle surface 
measured inwards, the stresses (i) to (iv) may be combined into the complete 
expression for the axial stress 


za|-7. + £7. - 7 +7] Tads 


ji 
+ - er for. — T,) sin @. ds — or D fcr, — T,) coso. as| (5) 


When the shape of the cylinder, and the temperature distributions, are 
given, it is a straightforward process to evaluate the coefficients in Equations 
(4), analytically or graphically, and to solve for X,, Y, and M,, and from these 
compute the corresponding stress. The axial stress can be found directly 
from Expression (5). 


Simple formulas can be derived for circular, elliptical, rectangular or other 
regular sections. The further development, with the exception of the next 
paragraph, will be confined to the circular case. 


Cylinder of Any Shape, but of Uniform Thickness, with Uniform but Unequal 
Inside and Outside Temperatures 


9. It is evident that since the temperature difference is constant around 
the cylinder, the curvature changes due to it in conjunction with the auxiliary 
end forces and couples will also be constant around the cylinder. Such a 
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curvature change can be annulled by a moment distribution M,. Clearly this 
is what is required to close the gap. It is easily found that 


M, = “Di eer =< 2). 


The extreme circumferential stress (6M,/h?), then, has the values* 
Ea(T, fab T;) 
~ 21 —o@) 
The complete axial stress is given by (ii) of Art 8, and has the extreme values 
+ $Ea(T. — T;) 
except near the ends, which, of course, are free of traction. It is noteworthy 
that these values are independent of the size, shape and thickness of the 
cylinder, depending only on the material and 
the temperatures occurring within it. 
It is known that these formulas hold for 
the special case of the circular tube (8, 
p. 373). 


Ve 





The Circular Cylinder 
10. When the cylinder has a mean radius 
a, the expressions on the left of the Equations 


(4) reduce to math(2M, + 3aX,), wanY,, ** 7 aa = 
2ma(M, + aX.) respectively, with axes as in M ¥{ M, 

Fig. 5; in the integrals on the right, ds is 7? ae 

replaced by adg, x by a sin ¢, and y by Fic. 5. 


a(l — cos@). The solution of the equations is 


xX = = (1+ a)a | fr. cos ddd — ef. — T;) cos 66] , 
D 27 26 
¥ = oat (1 ae a)a | |r. sin odd a ef. — T;) sin | ’ (6) 


M, = . (1+ eal — [Tacos sas + gf T-T0U + 2 cos | ‘ 


These formulas will give the bending moment, M, + X.y — Y.x, per unit 
axial length, explicitly in terms of the given temperature distributions. The 
direct stress due to X, and Y, will be small, of the order h/a, in comparison 
with the bending stress. 

When the temperature difference 7, — T; is as great as or greater than the 
mean temperature 7, = 3(7. + T7;), the contributions from the latter will 
be negligible, being of the order h/a, compared with the contributions arising 
from the temperature difference. 


* Positive signs pertain to the inside. 
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11. The integrals occurring in Expression (6) may be replaced by coefficients 
in the Fourier expansion of JT, and 7;. Let 4 
T. = Ayo + Aicosd® + Azcos2p + --- - @ 
+ Bysing + Besin2d +---- ie 
T; = Aj+ Aicosd + Aicos2¢ +---- 
+ Bising + Bysin2d +---- 

















Then 
27 27 27 
2wAy = [r. TA; = [.m 640 7B, ~ [tsi odo 
0 0 0 


and Aj, A{, Bi are similarly expressed in terms of T;. 
Thus X,, Y, and M,, and the stresses corresponding to them, are independent 
of the terms in cos 2q, sin 2q¢, and all higher harmonics of the temperatures. 


3 1 ; ‘ 
Replacing the factors C+ Sys Equations (6) can be written, 


kr 2 ahs 
D 
Ae os — a (lt + a)a(Ar — Al), 
D 
Y, = — tt + aa — Bi), (7) 


M. = 2 (1 + eal, — Af + A — Ad) . 


To obtain the corresponding extreme fibre stress, neglecting the direct stress 
due to X, and Y,, the moment M = M, + X.y — Y.x is multiplied by 6/h?. 
The result is 
eg age a aie B, — BY si } 8 
> 3 =~) t° 0 Ay 1) cos 6 + (B, 1) sin (8) 
whence the maximum follows easily. This is independent of the radius and 
thickness of the cylinder. 


The complete expression for the axial stress at the middle part of the length is 


Ea{-Tn + § (Ts T) + Ho + 48) 


| es 
+(A; = Al) 755 608 + (B, = Bt) Tq sin o} . (9) 
4 The last two terms will be relatively small. When they 


are neglected this stress becomes independent of a and h. 


12. As an example of thermal stress in a e 
circular cylinder, let the internal temperature ts 
be uniform, and let the external temperature 
exceed the internal by 37(1 + cos @). Then 
T is the maximum temperature difference, and 

X the distribution of temperature difference 
round the circumference is as indicated by 
the polar diagram, Fig. 6. It corresponds 

LT (1 + Cos d ) roughly to a tube heated from below, and 
containing a well mixed fluid. 
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Ay and A, are each equal to 37, and all other Fourier coefficients are zero. 
Then 


D Vf D 
X= —F(ltoas, Y, = 0,and MJ => (1 + o)alr . 


The moment M is given by 


> + ajar (1 - x). 


The maximum is at y = 0, and corresponds to an extreme fibre stress 
+EaT/2(1 — o@). 


To obtain a numerical result, take the internal temperature as 100° C., 
the maximum external temperature as 200° C., E as 3 X 107 Ib. per sq. in., 
a as 11.6 X 10-* per °C., o as 0.3 (as for steel). Then EaT/2(1 — @) is 
nearly equal to 25,000 Ib. per sq. in. 


The axial stress can be found from Expression (9). With the assumed 
distribution, there is a compressive stress in the hottest part (6 = y = 0) 
nearly equal to }EaT, or 27,000 lb. per sq. in. for the values taken. This 
is- the value for the stress at the inside. At the outside, the axial stress is 
9,000 Ib. per sq. in. tension. 


13. As a final example, let the outside of the cylinder be heated to a tem- 
perature 7, kept constant over a small arc of circumferential length A, and 
let the rest of the outside surface* and the whole inside surface be kept at 
zero temperature. This is a rough representation of a distribution that 
might occur in a welding process. 


It is easily shown that the maximum circumferential stress is given by 
3 HaF hk 


4i-¢ a 
Since this is proportional to \/a, it shows that localized heating does not 
produce large circumferential stresses. This is true, of course, only while 
there is no axial variation of temperature. 


As to the axial stress under such conditions, the parts (i) and (ii) of Art. 8 
are confined to the heated strip and occur throughout its length: (iii) and (iv) 
are distributions due to forces and moments localized on the ends of the heated 
strip. In the interior these will set up only very small stresses. In the 
interior, therefore, there will be an axial stress, given by—3EaT(1 + 2¢/h), 
confined to the heated strip, or a compressive stress varying from zero at the 
inside to EaT at the outside. 


The stress at the ends of the cylinder presents a separate problem, and can 
be dealt with approximately by methods given in another paper (2). It may 


* This disregard of the restriction to smooth temperature variations, referred to in Art. 3, is of 
no more importance here than its equivalent in plane stress systems with discontinuous boundary 
stress. The discontinuity conveniently replaces a variation as abrupt as is allowed, without 
seriously affecting the conclusions. 
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be shown that there would be local tensile circumferential stresses at the ends 
with a maximum value 2EaT/(3 + @). 

14. The thermal stress in a thick cylinder with a steady-state temperature 
distribution due to surface temperatures arbitrarily distributed round the 
inside and outside surfaces, but independent of the axial co-ordinate, can be 
reduced, as was shown by N. Muscheli§vili (5) and M.A. Biot (1), to an axial 
stress, and a dislocational stress corresponding to the stress specified above 
by means of X,, Y, and M,. In order to find the appropriate magnitudes of 
the dislocations it is necessary first to obtain the solution of the potential 
problem of finding the temperature distribution from the surface temperatures. 

By restricting the investigation described in this paper to the thin cylinder, 
the potential problem is avoided, and it becomes possible to find the stress 
and strain in terms of arbitrarily given surface temperatures for a cylinder 
of any shape and variable thickness. It should be said, however, that the 
potential problem presents no difficulties in the particular case of the uniform, 
thick circular cylinder. Explicit formulas for this case are given in another 


paper (3). 
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THE SPECIFIC HEAT OF COPPER FROM 30° TO 200° K.' 


By S. M. DocKErRTyY? 


Abstract 


An all-metal adiabatic vacuum calorimeter was used to determine the specific 
heat of copper over small temperature intervals from 30° to 200° K. The 
accuracy is considered to be within 0.05% for the greater part of the range. A 
curve is given showing the variation with temperature of the characteristic 
temperature, 0, for copper. 


Introduction 

Bronson, Chisholm, and Dockerty (2, 4) described the development of 
an all-metal adiabatic calorimeter and its use in precision measurements of 
the specific heat of copper from —80° to 100°C. The present paper deals 
with the construction of a similar type of calorimeter for much lower tem- 
peratures, used to determine the specific heat of copper from 30° to 200° K. 

The method, that of electrical heating under adiabatic conditions, has 
also been described (2,4). The general features can be seen from the diagram 
of the apparatus in Fig. 1. The calorimeter is cooled down to the lowest 
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Fic. 1. Adiabatic calorimeter. 
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point of the range at which one wishes to work, by immersing the container in 
a low temperature bath, the rate of cooling being speeded up somewhat by 
having hydrogen gas in the container. 


When the desired temperature is reached the container is evacuated. The 
specimen is heated so as to change its temperature in steps of 5° to 8°, the 
times of heating varying from 15 to25 min. Adiabatic conditions are obtained 
by keeping the temperature of the jacket as nearly as possible equal to that 
of the specimen, as indicated by the differential thermocouples. The tem- 
peratures are so adjusted, before and after the heating period, that the dif- 
ferences between the jacket and the specimen are less than 0.001° C. Thus 
the temperature of the jacket, as measured by the resistance thermometers, 
will be also the temperature of the specimen. The electrical energy supplied 
was measured, to within 0.02%, with a potentiometer and a standard cell 
that was recently calibrated in the laboratories of the National Research 
Council of Canada. The mass of the specimen was known to within 0.01%. 
Temperature differences were measured to within 0.02% for the greater part 
of the range, the precision being somewhat less than this at low temperatures. 


Jacket Apparatus 

Some difficulty was experienced with temperature gradients in the jacket of 
the older apparatus owing to the faulty distribution of heat loss. This 
was largely avoided in the present apparatus by using pure copper for the 
jacket, and by high evacuation which increased the thermal resistance between 
the jacket and the container about five times. Temperature gradients that 
were due to non-uniform heating were eliminated as before by winding separate 
heating coils on the sides and ends of the jacket; because of the high vacuum 
it was necessary to cement these coils in deep grooves. 


Differential thermocouples between the sides and ends of the jacket, as 
shown in Fig. 1, served to indicate temperature inequalities over the inner 
surface. These were put in from the outside and were left in permanently, 
so that adjustments could be made during the actual experiments. 


Two jackets were made by machining a 3-in. bar of copper. The first, 
which was used for most of the measurements above liquid air temperature, 
was about one and one-half times as large as that shown in Fig. 1, and dif- 
fered from it in that it was divided horizontally at the centre into two sym- 
metrical sections. This feature proved somewhat troublesome because it 
was found difficult to keep both sections at the same temperature. This 
equality was necessary because the resistance thermometers would indicate 
the temperature of the section in which their coils were situated. In addition, 
it was found that, because of the high vacuum, the size could be reduced con- 
siderably without affecting the accuracy of the results; this feature was 
desirable for the cooling to liquid hydrogen temperatures. A new jacket 
was therefore made as shown in Fig. 1. Care was taken to keep the end 
plug at the same temperature as the rest of the jacket, but a slight difference 
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would have no serious effect. Small leads of copper wire were used to prevent 
the production of cold spots on the jackets at low temperatures where the 
thermal conductivity of copper becomes quite high. 


Copper Specimen 

This consisted of a bar of commercially pure cold-rolled copper, a copper 
plug, and an 80-ohm manganin heating coil fitted together as shown in Fig. 1. 
The weights of the specimens used with the large and small jackets were 
about 650 and 320 gm. respectively. The amount of foreign material such 
as silk insulation, silk suspension, cement, and shellac was reduced to less 
than 0.1% of the total copper equivalent. 


Because of the high vacuum it was necessary to bind the heating coil to 
the plug with a thin coat of shellac. Small copper leads were used, and care 
was taken that these should be in good thermal contact with both the specimen 
and the jacket. 


The two-junction copper-constantan thermocouples were cemented in small 
holes in the specimen and the jacket respectively as shown. With a sensitive 
galvanometer these thermocouples would indicate temperature differences 
between the specimen and the jacket of less than 0.001° C. 


Resistance Thermometers 

Temperatures were measured with metal-stem resistance thermometers of 
platinum and lead placed in holes bored lengthwise in the jacket. To prevent 
heating and lagging effects it was necessary that the thermometers should not 
be evacuated with the rest of the system. They were therefore filled with 
helium gas at a pressure of one atmosphere, and sealed. A very satisfactory 
type of seal, which holds at low temperatures, is shown in Fig. 2. 


\ 
Vif 


Resistances were measured with a carefully cali- 
brated Miiller bridge maintained at a _ constant 
temperature in a thermostatic box. Small copper 
wire leads were used and inequalities in compensation 
were balanced out by a method involving two read- 


SOFT GLASS 
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TUBING 








ings on the bridge (7). 7 
The platinum coil was made of 0.002 in. wire eee 
(ro = 30 ohms) obtained from Johnson-Matthey, with 
a specified value for 6 (Equation (1) ) of 1.49. The teADS. 
lead coil was of 0.003 in. extruded wire (7) = 9 ohms) 5 
w 


obtained from the Baker Chemical Company, and 
was rated of high purity. 
RESISTANCE 


The platinum thermometer was calibrated in ice, 7 


steam, liquid oxygen, liquid hydrogen, and liquid 
helium in the usual way. The last two points 
proved to be of no interest from the point of view of 
calibration. In addition, a careful calibration was es 

made between — 183° C. and —210° C. by means _ Resistance thermometer. 


LEAD 
CAPILLARY 
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of an oxygen vapor pressure thermometer. Pure gas was obtained for this 
purpose by heating potassium permanganate in a vacuum. 
The values of the resistance ratio { R = - obtained for the various tem- 
0 
peratures are given in Table I. Repeated calibrations at some of these points 
showed no noticeable change over a period of a year. 
TABLE I 
RESISTANCE RATIO FOR PLATINUM THERMOMETER 











Temp., °C. R Temp., °C. R 
100.00 1.38976 —191.01 0.21407 
0.00 1.00000 —195.73 0.19377 
—182.98 0.24891 —210.41 0.13152 
—187.76 0.22823 —252.78 0.00957 
—188.66 0.22434 — 269.0 0.00438 


The lead thermometer was calibrated in ice and in liquid hydrogen, and 
indirectly for the range —180° to —210° C. by comparison with the platinum 
thermometer during the actual experiments. A slight zero shift was observed 
after the thermometer had been immersed in liquid air or in liquid hydrogen. 
The values of R obtained for the oxygen and hydrogen points were 0.29725 
at —182.98° C. and 0.03615 at —252.78° C. 

Resistance of Platinum. The international Temperature Scale (3) is defined 
in terms of the platinum resistance thermometer by the equation, 


a eS - 3 
Pie 8 (a5 aa) +71 ~ TF , (1) 


for 0° > ¢ > — 193°C. However, it has been found by Heuse and Otto 
(10), and by Keesom and Dammers (11), that temperatures derived by 
means of Equation (1) depart from the thermodynamic scale by as much 
as 0.05° C. at —80° C. and —140° C., and in opposite directions. 


Henning (9) has deduced that the resistances of different specimens of 
platinum should satisfy the relation, 

AR = a(R — 1) + 0(R —- 1)? , (2) 
where R is the resistance ratio 7/79... Onnes (13) has shown that this relation 
holds to — 215° C. with an accuracy in the derived temperatures of 0.02° C. 
Below this temperature no satisfactory relation has been found for comparing 
the resistances of different samples. 

Temperatures were calculated directly by means of equation (1) from the 
resistance measurements of the platinum thermometer, and also by comparison 
with a detailed standard calibration by means of Equation (2). The standards 
used were Henning’s thermometer No. 29 to —195° C., and Onnes’ No. 23’ 
from —195° to —210° C., the values of R for the standards being taken from 
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Leiden Comm. Supp. 58. Temperatures calculated by the two methods 
differed at some points by as much as 0.06° C. This gave rise to differences 
of as much as 0.1% in the temperature intervals. This is in agreement with 
the results of Heuse and Otto, and Keesom and Dammers. The method of 
comparison by Equation (2) was found to give close agreement with the 
calibration by means of the oxygen vapor pressure thermometer. The latter 
method was therefore adopted for the calculation of temperatures to — 210° C. 

Resistance of Lead. It has generally been found that the resistance ratios 
R, R’, of different specimens of lead can be compared by means of the relation, 


wm R-a ’ 
1—a 
which reduces to 
AR = a(R — 1). (3) 


Temperatures were first calculated by means of Equation (3), using as a 
standard a specimen calibrated in detail by Onnes (5) and deriving the 
constant a from the resistance at the hydrogen point. This method gave 
only approximate agreement with the calibration points at higher tem- 
peratures, and consequently the same Equation (2) that was used for platinum 
was also adopted for lead, with the constants a and b determined from the 
hydrogen and oxygen points. A further small correction term was added to 
allow for a discrepancy of 0.04° at —210°C. 


Experimental Measurements 

To test the reliability of the calorimeter under widely different experimental 
conditions some determinations of the specific heat, with the copper specimen 
above room temperature, were made; (A) with the calorimeter in a water 
bath, and (B) in a bath of liquid air. With the large calorimeter the values 
obtained from (B) were 0.15% higher than those from (A). In view of a 
similar correction in the previous investigation (4) it was thought that, with 
a liquid air bath, the correction would be proportional to the heat loss from 
the jacket. The measurements at temperatures above that of liquid air were 
made with this calorimeter. However, it was feared that at liquid hydrogen 
temperatures, where the heat capacity of copper is low, the above correction 
might have a larger value, and therefore the smaller calorimeter was con- 
structed for use in the lower temperature region. Results obtained above 
room temperature with this new calorimeter for Cases (A ) and (B) agreed 
within 0.02% of each other and within 0.05% of those obtained with the 
large calorimeter for Case (A). It was found that, when the differential 
thermocouples were kept balanced, the temperature of the specimen did not 
change by more than 0.001° C. in one hour. 

Measurements were made with the small calorimeter from 30° to 90° K. 
and two determinations were made at 93° K. and 133° K. At both these tem- 
peratures the results were 0.15% lower than those obtained with the large 
calorimeter. This showed that the correction factor did not vary as first 
supposed but was constant to 90° K. 
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Table II gives the values obtained for the specific heat of copper from 30° to 
200° K. Temperatures are given in the Kelvin scale (0° C. = 273.15° K.). 
The value of C, is given in joules per gram and in calories per mole (1 Cal..5° 
= 4.1835 electrical joules). The values of C, are calculated from the relation 

Cp— Cy = 1.60 X 10°C,2T . 
The values of the characteristic temperature, 6, were found from a table of 
Debye functions given by Beattie (1). A value of 5.9613 calories per mole 
was used for 3R. 

Since some of the temperature intervals were as great as 8° C., calculations 
were made to find the difference between the mean specific heat over the range 
in question and the true specific heat at the mean temperature. For all 
except the two lowest points the difference proved to be within the limits of 
experimental error. For the first and second points however the correction 
amounted to about 1%. This correction was quite accurately determined 


by a graphical method. 
TABLE II 


THE SPECIFIC HEAT OF COPPER 











Temp., C> Cy Co =5C, , 9 
nik. Joules per gm. | Cal. per mole | Cal. per mole | Cal. per mole 

28.64 0.02280 0.3475 0.0002 0.3473 314.0 
35.93 0.04440 0.6745 0.0005 0.6740 310.5 
42.58 0.06985 1.047 0.001 1.046 309.1 
50.13 0.09855 1.497 0.002 1.495 309.5 
59.24 0.1339 2.035 0.004 2.031 310.7 
67.21 0.1620 2.462 0.007 2.455 312.7 
74.64 0.1870 2.840 0.009 2.831 312.9 
87.45 0.2239 3.402 0.016 3.386 313.0 
87.88 0.2249 3.417 0.016 3.401 $13.2 
92.79 0.2365 3.594 0.019 3.515 313.8 
93.18 0.2373 3.606 0.019 3.587 314.0 
97.41 0.2463 3.743 0.022 3.721 314.9 
103.08 0.2579 3.919 0.025 3.894 S151 
108.51 0.2681 4.074 0.029 4.045 315.3 
193-73 0.2772 4.212 0.032 4.180 314.9 
119.38 0.2861 4.347 0.036 4.311 314.9 
125.42 0.2947 4.467 0.040 4.438 314.6 
131.29 0.3025 4.597 0.044 4.553 313.9 
132.97 0.3043 4.624 0.045 4.579 314.4 
137.48 0.3092 4.699 0.049 4.650 314.8 
144.24 0.3171 4.818 0.054 4.764 312.59 
151.08 0.3236 4.917 0.058 4.859 $12.3 
157.79 0.3293 5.004 0.063 4.941 312.0 
164.71 0.3345 5.083 0.068 5.015 312.0 
171.83 0.3399 5.165 0.073 5.092 310.4 
178.84 0.3443 S2a2 0.078 5.154 309.9 
186.33 0.3486 5.297 0.084 57243 309.6 
194.29 0.3528 5.361 0.089 S202 308.6 


Discussion of Results 
The accuracy of the results is limited in general by that of the temperature 
measurements. Above 80° K. it is considered to be within 0.05%. Below 
this temperature the accuracy is somewhat reduced because of uncertainties in 
the use of the lead thermometer. An error of more than 0.01° in a 5° interval 
in this range is considered unlikely. 
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Fig. 3 gives the values of @ plotted against the temperature: The corre- 
sponding values for the specific heat can be found from Beattie’s table. Above 
90° K. only one point lies off the curve by an amount corresponding to more 
than 0.05% in C,. The rise in the curve at 70° K. corresponds to about 
0.3% in C, and is thought to be well within the limits of error. The existence 


X- KEESOM ano ONNES 
O- PRESENT INVESTIGATION 


Characteristic temperature, 0 °K. 





o 50 100 /50 200 
Temperature, °K. 


Fic. 3. Characteristic temperature, 0, for copper. 


of a minimum at about 50° K. is in fair agreement with the work of Keesom 
and Onnes (12). Their points show a maximum deviation from the present 
curve of about 3%, which was the accuracy they had claimed. Both show 
an increase in the value of @ at lower temperatures. 

Above 110° K. the curve falls quite rapidly. In the previous work above 
200° K. it was found that the specific heat curve could be fitted with a constant 
value of 6 if a term, linear with the temperature, were added to Debye’s 
equation. This equation was found to hold, in the present case, to 100° K. 
with an accuracy of 1%. Below this temperature however the equation 
failed, since it amounted roughly to a linear correction in 6. 

The results agree with those of Griffiths and Griffiths (8) at room tem- 
perature but fall below their curve by 4% at 140° K. The results agree 
however at this temperature with the results of Eucken and Werth (6). 

In the previous work (2, 4) differences of 0.2% were found in the specific 
heat of specimens of specific gravities 8.8 and 8.9. The specific gravity of 
the specimen used in the previous work from —78° to 0° C. was 8.90, while 
that of the specimen used in the present work was 8.91. The two should 
have specific heats very nearly equal. Later unpublished work has shown 
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that the values previously obtained (4) are too low by 0.2%. When this 
correction is made the two agree to within 0.05%, both at room temperature 


and at 200° K. 
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